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ABSTRACT 

Neutron stars inspiralling into a stellar envelope can accrete at rates vastly exceeding the Eddington 
limit if the flow develops pressures high enough to allow neutrinos to radiate the released gravitational 
energy. It has been suggested that this hypercritical mode of accretion leads inevitably to the formation 
of stellar mass black holes during common envelope evolution. We study the hydrodynamics of this flow 
at large radii {R ^ Rns), and show that for low Mach number flows, in two dimensions, modest density 
gradients in the stellar envelope suffice to produce a hot, advection dominated accretion disk around 
the accreting object. The formation of outflows from such a disk is highly probable, and we discuss the 
impact of the resultant mass loss and feedback of energy into the envelope for the survival of the neutron 
star. Unless outflows are weaker than those inferred for well observed accreting systems, we argue that 
in most cases insufficient accretion occurs to force collapse to a black hole before the envelope has been 
ejected. This conclusions is of interest for black hole formation in general, for some models of gamma 
ray bursts, and for predictions of the event rate in future LIGO observations. 

Subject headings: accretion, accretion disks — hydrodynamics — binaries: close — stars: neutron — 
black hole physics — stars: evolution 



1. INTRODUCTION 

Common envelope evolution, in which the components 
of a binary system are engulfed by a common gaseous en- 
velope, is a brief but crucial phase in the formation of 
many compact binary systems. Drag forces, due to ve- 
locity differences between the orbiting components of the 
binary and the surrounding gas, work to shrink the binary 
orbit, while the potential energy released acts to expel the 
common envelope. Extensive numerical work (e.g. Boden- 
heimer & Taam 1984; Livio & Soker 1988; Taam, Boden- 
heimer & Rozyczka 1994; Sandquist et al. 1998) has only 
partially succeeded in reducing the large uncertainties sur- 
rounding the efficiency of this process (Iben & Livio 1993; 
Livio 1994; Rasio & Livio 1996). 

For neutron stars, an equally fundamental uncertainty 
concerns the accretion rate during the common envelope 
phase. At the densities of a typical giant envelope, the 
Bondi-Hoyle accretion rate A/bh for an inspiralling neu- 
tron star would be extremely large, typically of the order 
of '^ 1 M0yr~^, and often much larger still. This is many 
orders of magnitude in excess of the Eddington rate^ of 
^ 10^* Af0yr^^ — obtained by equating the outward ra- 
diation pressure with gravity — leading to starkly differ- 
ent predicted outcomes. If M is limited to the Eddington 
rate (or is even several orders of magnitude higher) , then 
ejection of the stellar envelope will occur when the total 
accreted mass is Mace ^ Mq, and the neutron star will 
survive. Conversely, uninhibited accretion at close to A/bh 
will lead almost inevitably to collapse to a black hole. 

Recent work has tended to suggest that A/bh is the ap- 
propriate rate for neutron star accretion in a stellar enve- 



lope (see e.g. Bethe & Brown 1998 and references therein). 
At sufficiently high accretion rates, the Eddington limit 
becomes irrelevant because photons are trapped and ad- 
vected inward with the flow (Rees 1978; Begelman 1979; 
Blondin 1986). A hot envelope develops around the neu- 
tron star, and eventually high enough densities and tem- 
peratures are reached that the accretion energy can be 
radiated away by neutrinos. This hypercritical mode of 
accretion has been considered in the context of fallback 
onto neutron stars in supernovae (e.g. Colgate 1971; Zel- 
dovich, Ivanova & Nadczhin 1971; Chevalier 1989; Houck 
& Chevalier 1991), as a component of models for gamma- 
ray bursts (Popham, Woosley & Fryer 1999), and applied 
to the problem of accretion in a stellar envelope (Chevalier 
1993, 1996; Brown 1995; Fryer, Benz & Herant 1996). A 
striking consequence of ubiquitous hypercritical accretion 
would be that neutron star-neutron star binaries cannot 
form from massive binaries via the usual common envelope 
evolution route, which would instead lead to black hole- 
neutron star binaries. The observed population of neutron 
star-neutron star binaries would have to form instead via a 
rarer channel in which the binary progenitors are of nearly 
equal mass (Brown 1995), and would be an order of magni- 
tude rarer than black hole-neutron star binaries (Bethe & 
Brown 1998). The robustness of this conclusion is of great 
interest, since these binaries are both potential gamma-ray 
burst progenitors, and the most promising targets for the 
early LIGO detection of gravitational waves (Abramovici 
et al. 1992). 

The velocity of the neutron star relative to the com- 
mon envelope ensures that the accretion flow cannot be 
spherically symmetric, at least at large radii of the order 
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of the accretion radius. Moreover, gradients in the enve- 
lope structure across the accretion radius introduce angu- 
lar momentum into the flow, which modifles its properties 
even at small radii, via the formation of a rotationally sup- 
ported disk. In this paper we investigate how these com- 
plexities affect hypercritical accretion via hydrodynamic 
simulations of the outer parts of the accretion flow. Our 
approach thus complements the consideration of angular 
momentum presented by Chevalier (1996), though as dis- 
cussed later we draw different conclusions as to its likely 
importance. 

The plan of this paper is as follows. In §2 we briefly 
review the properties of spherical hypercritical accretion 
that are relevant for the common envelope application, 
and in §3 we outline the numerical methods used. The 
principal computational limitation is the restriction to two 
dimensional simulations. Results for flows with zero and 
non-zero angular momentum are presented in §4 and §5, 
and the implications and remaining uncertainties discussed 
in §6. 

2. HYPERCRITICAL ACCRETION 

The simplest approximation to the accretion flow in 
a common envelope is to describe it as Bondi-Hoylc- 
Lyttleton accretion from a uniform medium onto a point 
mass. The properties of such flows have been extensively 
studied, both analytically (Hoyle & Lyttleton 1939; Bondi 
& Hoyle 1944), and numerically (e.g. Ruffert 1997; Kley, 
Shankar & Burkert 1995; Benensohn, Lamb & Taam 1997, 
and references therein). For an accretor of mass M, the 
critical radius is the accretion radius^, 



Ra 



2GM 



(1) 



where Woo and Cqo are the velocity and sound speed far 
upstream of the accreting star. Material falling into a 
cylinder with this radius is accreted, so that. 



Mbh = ttRIpooV^ 



(2) 



where p^o is the density of the gas far upstream of the ac- 
cretor, and the last factor is included as it provides better 
results for low Mach number accretors. The timescale for 
setting up the flow is of the order of the sound crossing time 
at the accretion radius, ta = Ra/coo- Generally, since we 
are interested in inspiral through a pressure supported at- 
mosphere, we expect Mach numbers of order unity, though 
this will vary somewhat at different radii in the star. 

Numerical values for Ra and Mbh require specification 
of a stellar model. Fryer, Benz & Herant (1996) tabulate 
values for several giant and main sequence stars (the giants 
are of greater interest here, though somewhat analagous 
considerations apply to neutron star-main sequence colli- 
sions, for example in globular clusters). For example, for 
a 20 Mq giant Ra ^ 10^^ cm, while Mbh varies from 
10^ MQyr~^ on upwards. For a more distended 10 Mq 
model lower accretion rates of 10~^ MQyr~^ are of course 
expected in the outer regions. These accretion rates can 
only be realized in the hypercritical regime, since the Ed- 
dington limit on the luminosity. 
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where nip is the proton mass and ax the Thomson cross- 
section, corresponds to a few xlO^* MQ-yr^^. 

For hypercritical accretion to occur the timescale for 
photons to diffuse out of the flow must exceed the timescale 
for them to be advcctcd inwards. For a spherical flow 
whose opacity Kcs is dominated by electron scattering, the 
optical depth is 



p{R)KcsdR, 



(4) 



and the diffusion time is roughly. 
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For a flow in free-fall, where Vr — ~y^2GM/R, equating 
idiff to the free-fall time R/vr yields an estimate for the 
trapping radius (Blondin 1986), 



_ 1 Mc^ 2GM 
^*^^P-2LEdd c2 ■ 



(6) 



In addition to the photons being trapped, we must also 
be able to eventually lose the accretion energy from the 
system. For a black hole, this is possible at any accretion 
rate since the energy can be invisibly advected across the 
horizon (e.g. Narayan, Garcia & McClintock 1997). In a 
neutron star system, high accretion rates are required be- 
cause the pressure at the base of the accretion flow must 
be able to become large enough to allow neutrino losses 
(primarily from electron positron annihilation into neu- 
trino pairs) to balance the accretion energy, GMM/Rns- 
In the spherically symmetric case, this requirement usually 
requires the formation of a hot envelope at small radii, in 
which Vr is subsonic, bounded from free falling gas by a 
shock at radius i?sh- Within this inner region, the den- 
sity and pressure scale as a radiation dominated 7 = 4/3 
{n = 3) polytrope (e.g. Brown 1995; Fryer, Benz & Herant 
1996), 
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A detailed calculation (Houck & Chevalier 1991; Brown 
1995) flnds for R^h, 



i?sh '^ 2.6 X 10^ cm 



M 



-0.37 



Mq yr-1 



(8) 



ut 



(3) 



Requiring that i?sh < ^trap gives a limit on the accretion 
rate, of roughly Mhypcr ^ 10"'* Moyr"^. Even allow- 
ing for a generous margin of error in this crude estimate, 
for example due to a non-spherical geometry, this is sub- 
stantially less than the Bondi-Hoyle rate in an envelope, 
which should thus be safely in the hypercritical regime. Of 



course we also require in the common envelope case that 
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Rsh < Ra, but for spherical accretion and Ra ^ 10^^ cm 
this is an easily satisfied constraint, Af > 10"^ AfQyr"^. 

If angular momentum and other possible complications 
can be neglected, the scenario for the inspiral of a neutron 
star then has three stages. First, the neutron star accretes 
at less than the Eddington rate as it encounters the tenu- 
ous outer atmosphere of the companion. Negligible mass 
is accumulated in this phase. As the density grows, the 
implied Bondi-Hoyle rate first rises above AfEddj but is 
below Mhyper- In this regime outward diffusion of pho- 
tons is effective at limiting accretion, again to negligible 
levels. Finally, the Bondi-Hoyle rate rises well above the 
critical value for hypercritical accretion, and gas starts ac- 
creting freely onto the neutron star. The outcome then 
depends critically on how much mass can be accumulated 
before enough energy has been lost to unbind the envelope. 
Brown (1995) estimates this energy as. 



E^3<vL> AM 



(9) 



where AM is the accreted mass, and < w^ > is the av- 
erage mean velocity squared during the hypercritical ac- 
cretion phase. For typical binding energies and velocities, 
AAf ^ 1 Mq, sometimes vastly so, and the mass accreted 
will probably exceed the maximum mass for a neutron star 
and force collapse to a black hole. 

3. NUMERICAL METHOD 

We investigate how angular momentum could affect the 
above scenario of hypercritical accretion by studying the 
inviscid, purely hydrodynamic behavior of infalling gas at 
large radii, R ^ Rns- The hypercritical regime corre- 
sponds to the complete dominance of advection of photons 
over diffusion, so that the effects of radiation transport are 
negligible. We also neglect magnetic fields, which is po- 
tentially less forgivable. To attain the required dynamic 
range in three dimensional calculations including magnetic 
fields unfortunately remain difficult. 

3.1. 2D simulations 

We have investigated two dimensional simulations both 
in spherical polar geometry (with v parallel to 6 = and 
axisymmetry in the (/) direction), and in cylindrical geome- 
try (z, R, (jj). The spherical polar simulations represent the 
most faithful two dimensional representation of the flow, 
but necessarily exclude the possibility of studying angular 
momentum or density gradients in the ambient medium. 
For such zero angular momentum flows, we found only that 
the analytic Bondi-Hoyle estimate of equation (0) provides 
a good estimate of the accretion rate (for low Mach num- 
bers the actual accretion rate is somewhat higher than the 
estimate, but the difference is not so great as to qualita- 
tively affect the result). 

For these reasons (and following Benensohn, Lamb & 
Taam 1997 and most other two dimensional calculations) 
we focus on the cylindrical polar calculations. The neutron 
star is represented as a point mass at i? = 0, surrounded 
by a fixed computational mesh in {R, 0). We use uniform 
zoning in the (f> direction, and choose the radial grid such 
that Ri+i = (3Ri with /3 > 1 a constant. This amounts 
to choosing cells that have the same shape, in our case 
roughly square, at all radii. 



Cylindrical polar simulations in which z is the ignored 
co-ordinate allow for angular momentum, but correspond 
to the rather unphysical assumption that the central mass 
moves through a thin sheet of gas with zero gradients in 
the perpendicular direction. This may lead to qualita- 
tive errors, for example in the strength and prominence 
of 'flip-flop' instabilities at high Mach numbers (compare 
e.g. Livio 1992; Benensohn, Lamb & Taam (1997); Ruf- 
fert (1997)), though untangling the effects of geometry and 
resolution in the various calculations is difficult. We will 
not be interested here in transient disks or instabilities in 
the wake, where these issues are most worrisome. 

We assume that the gas is dominated by radiation pres- 
sure out to the outer boundary of the simulation at a few 
Ra- The equation of state can then be modelled using a 
simple adiabatic relation, 



P= h - l)e, 



(10) 



where e is the energy per unit volume and we take 7 = 
4/3, corresponding to a radiation dominated gas where 
p=l/3e. 

The calculations use the ZEUS-3D code developed by 
the Laboratory for Computational Astrophysics (Clarke, 
Norman & Fiedler 1994). ZEUS is an Eulerian finite 
difference code which employs an artificial viscosity to 
handle shocks. The algorithms and design of the code 
are closely similar to those detailed by Stone & Norman 
(1992a, 1992b). 

3.2. Boundary conditions 

For the outer boundary condition, we impose inflow of 
gas at velocity Wqo and sound speed Coo at the outer bound- 
ary. Rout, for —tt/2 < (p < 7r/2. Over the remainder of the 
outer boundary outflow boundary conditions are specified, 
implemented as simple continuation of fluid variables on 
the grid into the boundary zones. We take a Mach number 
A4 — 1.5, appropriate to the generally mildly supersonic 
accretion flows in common envelopes. 

An inspiralling neutron star encounters radial gradients 
in poc, Coo and ?;oo across the accretion radius. In most 
cases the density gradient is the principal effect (e.g. Fryer, 
Benz & Herant 1996), which we model as a simple expo- 
nential, 

Poo oc e^''^'-/^" (11) 

where Ar is the radial distance in the common envelope of 
the unperturbed medium from the neutron star orbit, and 
ep measures the strength of the gradient across the accre- 
tion radius. Values of ep for various stars vary substan- 
tially, and have been tabulated by Fryer, Benz & Herant 
(1996). 

The inner boundary condition is reflecting (vr = at 
R = -Rin), allowing for the formation of a pressure sup- 
ported inner envelope for zero angular momentum accre- 
tion. This implicitly assumes low accretion rates where 
the shock radius Rgh is large enough to exceed -Rin. For 
runs with angular momentum, the choice of inner bound- 
ary condition is less important as angular momentum pro- 
vides substantial support against gravity at small radii. 
Our boundary condition then amounts to assuming that 
radial flow through the disk is slow - the plausibility of 
this can be verified post facto by studying the properties 
of the disks formed in the calculation. 
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4. DISK FORMATION 

Figure 1 shows results for a series of Bondi-Hoyle ac- 
cretion simulations in which the density gradient in the 
ambient medium was varied from ep = to ep = 0.4. All 
the calculations used a grid with n^ — 144 and ur — 160. 
The inner boundary was at i?in = -Ra/60 and the outer 
boundary at i?out = 4i?a, giving a grid with /3 ~ 1.03. 
The calculations were run until t = 32ia, where the time 
unit ta is the sound crossing time of the accretion radius 
Ra- We plot in Fig. 1 only the inner region of the accretion 
flow. The ep ~ 0.2 run was also recomputed at modestly 
higher resolution (n^ = n^ = 200) until t — lOOfa, in or- 
der to ascertain how quickly mass accumulated in the outer 
parts of the disk at a later time. No qualitative changes 
were observed to occur during this longer simulation. 

In the absence of density gradients in the ambient gas, 
the structure of the flow resembles closely that seen in 
simulations of Bondi-Hoyle accretion with the same pa- 
rameters and an absorbing inner boundary condition. A 
pressure supported, roughly symmetric envelope has de- 
veloped around the central object, and this displaces the 
bow shock upstream into the flow. At this resolution and 
Mach number, the flow is found to be only rather weakly 
transient. Increasing the density gradient first leads to 
an asymmetric displacement of the bow shock, followed 
for more extreme density gradients by the formation of a 
clear disk in the inner regions of the simulation. For the 
parameters adopted here, Cp ^ 0.2 suffices to create a clear 
and persistent disk surrounding the central object. 

Figure 2 plots the velocity and density fields correspond- 
ing to the images in Fig. 1. Over this radial range there 
is a large variation in velocity, which we normalise to the 
local Keplerian value Vk — y/{GM/R). For the zero den- 
sity gradient case, the outer flow is clearly not cylindrically 
symmetric, but within ~ 0.1i?o of the central object the 
envelope is both reasonably symmetric and characterised 
by only small rotational motions. Conversely, a clear disk 
is produced in the calculations with stronger density gra- 
dients, ep > 0.2, the flow here is disk- like out to at least 

Ra/2. 

Figure 3 shows how the volume averaged rotational ve- 
locity from the simulations, normalised to Vk, varies as a 
function of radius. Negligible rotation is seen in the ep = 
run, with < v^ > /vk < 0.1 at all radii, while the runs 
with density gradients all show the signature of a disk in 
which both pressure gradients and rotational support are 
significant. For the two runs with the largest gradients, 
a sizeable disk extending out to almost Ra, with typical 
angular velocity fi w Vtk/2, is clearly produced. 

5. DISK EVOLUTION 

For hypercritical accretion, the disks formed around the 
neutron star will be advection dominated. The proper- 
ties of advection dominated flows have been extensively 
studied, both at the high accretion rates of interest here 
(e.g. Begelman & Meier 1982) and more recently at lower 
accretion rates (e.g. Narayan & Yi 1994, 1995). These 
disks are hot, geometrically thick, and only weakly bound 
to the accreting object. The presence of a disk may affect 
the outcome of common envelope evolution in two ways, 
via a modification of the accretion rate onto the neuton 
star or via feedback of energy into the stellar envelope. 



The most direct potential influence of a disk is via the 
accretion rate. As infall continues, the disk will reach a 
quasi-steady state in which the rate of infall onto the outer 
disk balances the rate of disk accretion. This rate need not 
be the Bondi-Hoyle rate, and in principle could be lower 
(for an illustration, in our inviscid simulations, where the 
viscosity is very low, we eventually reach a steady state 
where there is close to zero ongoing accumulation of mass 
in the disk). We estimate whether this is important be- 
low, and show that for the usually considered values of 
q;ss7 the Shakura-Sunyaev (1973) viscosity parameter, the 
disk is probably able to transport mass inwards at the 
Bondi-Hoyle rate. 

The presence of a disk also makes the formation of out- 
flows or jets probable (e.g. Livio 1999). Strong outflows 
could themselves affect the fate of the system by reducing 
the accretion rate onto the neutron star. Much weaker 
jets, if they arose from deep in the neutron star potential 
well, would still provide an important additional energy in- 
put into the envelope, and shorten the common envelope 
phase. 

5.1. Accretion rate 

The accretion rate through the disk can be estimated 
from the measured disk mass and the inferred viscous 
timescale. For a viscosity parameterized via the Shakura- 
Sunyeav (1973) form, v = assCg/f^fc, and at radius R the 
viscous timescale is. 



t,. 



i?2 



R ^k 

assCg 



(12) 



The mass in the disk will drain onto the central object on 
the viscous timescale at the outer radius, with an accretion 
rate Mdisk ~ Mdisk/t^- The value of ass appropriate to 
thick, radiation dominated disks is unknown. If MHD in- 
stabilites (Balbus & Hawley 1991) arc responsible for the 
viscosity, then simulations of thin gas pressure dominated 
disks find typically that ass ~ 10~^ (Stone et al. 1996; 
Brandenburg et al. 1996). There remain large uncertain- 
ties in the theoretical expectation for disks of the kind 
that we are interested in. However, taking ass = 10~^, 
and evaluating the accretion rate from the disk mass and 
run of sound speed obtained at t — lOOia in the long du- 
ration ep = 0.2 run, we find that Mdisk ^ Mbh- This is 
a crude estimate, which in particular ignores completely 
the expansion of the disk expected from viscous evolution 
(Lynden-BcU & Pringle 1974). As a result, the value Afen 
should be regarded as an upper limit to the mean accre- 
tion rate through the disk. Furthermore, at least at small 
radii the accretion is unlikely even to be steady (Chevalier 
1996), since the much slower fall-off of pressure with radius 
in an advection dominated disk as compared to a spheri- 
cal envelope (p oc R~^'^ as compared to p oc R^^) requires 
non-steady accretion to reach the extreme conditions re- 
quired for neutrino emission. However, it suggests that 
there is no strong reason to believe that, in the absence of 
outflows, the formation of a disk creates an insurmount- 
able bottleneck to rapid accretion. 

5.2. Outflows and jets 

Most accretion disk systems are observed to generate 
jets or less well collimated outflows. Although a thick disk 
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such as that formed in the simulations has a relatively 
small density contrast between the equatorial plane and 
the polar regions, advection dominated flows are also hot 
enough that the gas is only weakly bound to the central 
mass. Detailed solutions show that the Bernoulli constant, 



deposition after AMjet of mass has been ejected is just, 
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(13) 



which measures the energy the gas would possess if adi- 
abatically moved to infinity, is positive for an often wide 
range of angles close to the poles (Narayan & Yi 1994; 
1995). Although the outcome depends additionally on the 
outer boundary conditions and the detailed physics of out- 
flow generation, this positivity of Be is likely to imply that 
outflows are a generic feature of advection dominated disks 
( Narayan & Yi 1995; Blandford & Begelman 1999). For 
our purposes, we can distinguish two extreme possibili- 
ties, in which outflows are either self-similar or generated 
exclusively from the inner disk at i? ~ i?ns- 

5.2.1. Outflows 

Self-similar outflows from advection dominated flows 
represent the model considered by Blandford & Begelman 
(1999). In this case the fraction of accreting mass lost in 
the outflow is the same for each decade in disk radius, so 
that the remaining mass accretion rate through the disk 
decreases inwards as M oc i?" with < n < 1. The appro- 
priate value of n, which within this model measures the 
efficiency with which mass is ejected, has not been deter- 
mined for any realistic thick disk model, though it is has 
been suggested that n is large for simulations of convection 
in thick disks (Stone, Pringle & Begelman 1999). Empir- 
ically n must be n ~ 1 if this model is to be successful in 
explaining the extremely low radiative efficiencies of disks 
around black holes in low luminosity galactic nuclei (e.g. 
Reynolds et al. 1996). In the common envelope case, the 
disk extends over an extremely large range of radii, from 
the neutron star surface at i? ~ 10^ cm out to of the order 
of the accretion radius at i? ^ 10^^ cm. As a result, the 
consequence of outflows that are inefficient at removing 
angular momentum would be to greatly reduce the accre- 
tion rate onto the neutron star. Cygnus X-2 may be an 
example of a neutron star that has survived an episode of 
super-Eddington mass transfer (though here during ther- 
mal timescale mass transfer rather than common envelope 
evolution) without accreting a significant fraction of the 
transferred mass (King & Ritter 1999; King & Begelman 
1999). If this is indeed the case, it provides some support 
for the scenario of strong outflows under physical condi- 
tions analagous to those encountered during common en- 
velope inspiral. 

5.2.2. Jets 

Alternatively, outflows may arise predominantly from 
the inner disk. In this case, the energy feedback into the 
common envelope from outflows originating deep in the 
neutron star potential well could be highly significant. If, 
as observations of jet systems suggest (Livio 1999), the 
jet is launched with a velocity roughly equal to the lo- 
cal Keplerian velocity, vf^^. ~ GMns/Rns, then the energy 



E^ 



GM^sAMi 



'jet — Q^jot" 



jet 



Rn 



(14) 



where ajct is an efficiency factor that will depend on the 
specific jet model. Ejecting the common envelope requires 
an energy deposition of around 2 x 10^^ ergs (Brown 1995), 
which could be achieved for a mass loss in a jet as low as. 



AMjct ~ 5 X 10 



-6„-l 



"je; Mq. 



(15) 



The energetic feedback from such an outflow could thus be 
important for hastening the ejection of the envelope, even 
if the mass loss itself was far too small to significantly 
impact the accretion rate. Indeed for this type of jet a 
large energy deposition could only be avoided if either the 
accretion rate was ^ Mbh, or if the efficiency of the jet 
production or coupling to the envelope was extremely low. 
A low efficiency of coupling to the envelope could arise if 
the jet was extremely well coUimated and able to escape 
the star entirely, which remains a possibility. However, 
even for a collimation angle of ^ 10^^, which is typical 
for many jet sources, the ejection of ~ 0.2Mq is sufficient 
to unbind the envelope. If jets are able to deposit energy 
into the envelope then an immediate consequence would 
be that the outcome of common envelope evolution should 
depend on the depth of the potential well at the surface 
of the inspiralling compact object - more compact objects 
should lead to a higher value of the common envelope effi- 
ciency parameter acE and eject the envelope more easily. 

6. DISCUSSION 

In this paper we have considered the probable fate of 
neutron stars during common envelope evolution. Nu- 
merical simulations of the hydrodynamics of Bondi-Hoyle 
accretion at large radii from the neutron star show that 
modest density gradients, typical of those expected in gi- 
ant envelopes, lead to persistent rotationally supported 
disks around the accreting object. The presence of angular 
momentum means that spherically symmetric treatments 
of hypercritical accretion are likely to be a poor approxi- 
mation for studying the the common envelope phase. In 
particular, whether the neutron star accretes enough mass 
to force it to collapse to a black hole prior to loss of the 
envelope depends entirely on the subsequent evolution of a 
thick, advection dominated, accretion disk. Little is known 
about the behavior of such disks, opening plenty of room 
for legitimate dispute as to the outcome. 

Simple estimates suggest that the disks formed in the 
simulations are sufficiently hot and thick to support accre- 
tion rates which could be as large as the inferred Bondi- 
Hoyle rate in the common envelope phase, provided only 
that the disk viscosity is not surprisingly small. The pres- 
ence of a disk does not create a serious bottleneck in the 
accretion flow at _R '-^ Ra- However, more complex and 
uncertain physics is likely to come into play close to the 
neutron star. Jets are observed almost universally from 
accreting systems possessing thick, hot accretion disks, of 
the kind envisaged here. In the case of SS433, which ar- 
guably is the closest prototype to the physical condition 
of a neutron star in a common envelope, the mass out- 
flow in the jets is extremely strong (e.g. Watson et al. 
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1986). During common envelope evolution, a jet would 
allow for a strong feedback of energy from close to the 
neutron star surface into the stellar envelope, leading to a 
more rapid ejection of the envelope than would be possible 
from gas and gravitational drag alone. Since the fate of 
the neutron star depends on a sensitive balance between 
the rate at which it accretes and the rate at which energy 
is deposited into the common envelope (Brown 1995), this 
would improve the chances of neutron star survival both 
by lowering the accretion rate and by reducing the epoch 
of common envelope evolution. In general, jets must be 
avoided at all costs if the neutron star is to be able to ac- 
crete a large mass during inspiral. It is also possible that 
enough of the radiated neutrino energy could be absorbed 
at a larger radius to drive explosions, as discussed in the 
context of spherical accretion models by Fryer, Benz & 
Herant (1996). 

Observationally, several binary pulsars are known whose 
properties would be consistent with the neutron star 
having survived a phase of common envelope evolution. 
Camilo et al. (1996) identify four pulsars which have rela- 
tively large companion masses (in excess of 0.45 Mq), and 
which do not follow the eccentricity-orbital period relation 
expected for lower mass binary pulsars (Phinney 1992). 
These systems are likely to have undergone deep common 
envelope evolution (van den Heuvcl 1994), although prob- 



ably the companions were of rather lower mass (1-6 Mq) 
than we have been discussing here. Nonetheless, what is 
striking is that these neutron stars appear to have not 
only survived, but are now observed to be rotating with 
spin periods that are rapid, yet significantly slower than 
pulsars believed to have been spun-up via disk accretion. 
This is consistent with the scenario argued here, in which 
common envelope evolution involves accretion rates that 
are vastly super-Eddington yet still insufficient to force 
collapse of the neutron star to a black hole. It is also con- 
sistent with the suggestion of Brown, Lee & Bethe (1999), 
who point out that the inferred masses of black holes in 
transient sources (which appear to be a large fraction of 
the mass of the immediate progenitor) imply that hyper- 
critical accretion onto black holes in supernovae must be 
reasonably efficient. In supernovae, mass must be physi- 
cally ejected from the system to avoid being accreted even- 
tually, whereas in the common envelope case a modestly 
lowered accretion rate can suffice to enable the neutron 
star to survive until the envelope is lost. 
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Fig. 1. — Density (on a logarithmic scale) for different values of the density gradient €p. Top left ep = 0, top right ep = 0.1, lower left 
ep = 0.2, lower right ep = 0.4. The images measure 4Ra on a side. 
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Fig. 2. — Density and velocity fields in the inner region of the flow. The panels represent the same sequence of increasing ep as Fig. 1, with 
contours plotted at Alog(p) = 0.1. The length of the vectors has been scaled to the local Keplerian velocity Vf^ oc R~^'^. 
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Fig. 3. — Degree of rotational support for different values of £p. Plotted is the volume averaged rotational velocity -u^, normalised to the 
local Keplerian velocity v^, as a function of the radius R in units of the accretion radius -Ra- The solid line shows results for €p = 0, the 



dotted line i 



0.1, the short dashed line ep = 0.2, and the long dashed line €p = 0.4. 



